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Abstract. The quark-hadron duality is studied in a systematic way for polarized and unpolarized 
structure functions, by taking into account all the available data in the resonance region. In both 
cases, a precise perturbative QCD based analysis to the integrals of the structure functions in the 
resonance region has been done: non perturbative contributions have been disentangled, the higher 
twist contributions have been evaluated and compared with the ones extracted in the DIS region. A 
different behavior for the unpolarized and polarized structure functions at low Q 2 has been found. 
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The structure and the interaction of hadrons is generally described by two different but 
complementary approaches: the quark-gluon context at high energy, where the quarks 
are asymptotically free, and the meson and baryon description at low energy, where the 
effects of the confinement are large. In some specific cases where the natural description 
in terms of hadrons should be applied, the quark-gluon description can be successfully 
also used. This evidence is called quark-hadron duality and it was introduced by Bloom 
and Gilman who noted a relationship between the nucleon resonance region and the 
deep-inelastic- scattering (DIS) one yj]. 

A quantitative analysis of the Q 2 dependence of quark-hadron duality in both polar- 
ized and unpolarized ep scattering is presented. All current data in the resonance region, 
1 < W 2 < 4 GeV 2 , have been taken into account. For the unpolarized case it has been 
used the data obtained at Jefferson Lab in the range 0.3 < Q < 5 GeV |2fl, and the data 
from SLAC (J3J] and references therein) for Q 2 > 4 GeV 2 . For the polarized case there 
are only few experimental data in the resonance region. One set is part of the El 43 data 
1 4], and it corresponds to Q 2 = 0.5 and 1.2 GeV 2 . Another set is the one from HERMES 
1 10] in the range 1.2 < Q 2 < 12 GeV 2 . 

The full procedure of the analysis is described in [7]. The quark-hadron duality 
in DIS is studied by considering the ratio of the integrals of the structure functions 
integrated in a defined x-range, corresponding to the W range of the resonance region. 
The structure function in the numerator is evaluated using the experimental data in the 
resonance region, while the one at the denominator is calculated from parametrizations 
that reproduce the DIS behavior of the data at large Q 2 . The ratios T^pd and R po \ have 
been calculated in unpolarized and polarized cases, respectively. 

In order to understand the nature of the remaining Q 2 dependence that cannot be 
described by NLO pQCD evolution, the effect of target mass corrections and large x re- 
summation have been studied. The analysis was performed by using x as an integration 
variable, which avoids the ambiguities associated to the usage of other ad hoc kinemati- 



cal variables. Standard input parametrizations with initial scale Q 2 = 1 GeV 2 have been 
used. Once both effects have been subtracted from the data, and assuming the validity of 
the twist expansion, one can interpret any remaining discrepancy of the ratio from unity 
in terms of higher twist. 

The Target Mass Corrections (TMC) are necessary to take into account the finite 
mass of the initial nucleon. They are corrections to the leading twist (LT) part of the 
unpolarized structure function F2. 

Large x Resummation (LxR) effects arise formally from terms containing powers of 
ln(l — z), z being the longitudinal variable in the evolution equations, that are present 
in the Wilson coefficient functions C(z). The logarithmic terms in Cns(z) become very 
large at large x, and they need to be resummed to all orders in as [&]. 

All the effects are summarized in Fig[H where the ratio between the resonance 
region and the 'DIS' one is reported for the unpolarized and for the polarized case: 
the numerator is obtained from the experimental data, while the denominator includes 
the different components of the present analysis, one by one. 

For unpolarized scattering it has been found that TMC and LxR diminish considerably 
the space left for HT contributions. The contribution of TMC is large at the largest 
values of Q 2 because these correspond also to large x values. Moreover, the effect of 
TMC is larger than the one of LxR. Similarly, in polarized scattering the inclusion of 
TMC and LxR decreases the ratio R^ ol - However, in this case these effects are included 

almost completely within the error bars. Clearly, duality is strongly violated at Q 2 < 1.7 
GeV 2 . The present mismatch between the unpolarized and polarized low Q 2 behavior 
might indicate that factorization is broken differently for the two processes, and that the 
universality of quark descriptions no longer holds. The discrepancy from unity of the 
ratios already presented is interpreted in terms of HTs. In Figs. I2l3l the question of the 
size of the HT corrections is addressed explicitely. For F%, they are defined as: 

H(x, (?) = <£- (F 2 - 2 2 ) - FV) ; C„ w = J£f*> - Q^T^ZIlL 

^2 (x/Q z ) r 2 

(1) 

A similar expression is assumed for g\. Cht is the so-called factorized form obtained 
by assuming that the Q 2 dependences of the LT and of the HT parts are similar and 
therefore they cancel out in the ratio. Although the anomalous dimensions of the HT 
part could in principle be different, such a discrepancy has not been found so far in 
accurate analyses of DIS data. The HT coefficient, Cht has been evaluated for the 
three cases listed also in FigUl namely with respect to the NLO pQCD calculation, 
to NLO+TMC and to NLO+TMC+LxR. The values of 1 +C H t/Q 2 are plotted in FigEl 
(left panel) as a function of the average value of x for each spectrum. One can see that 
the NLO+TMC+LxR analysis yields very small values for Cht in the whole range of x. 
Furthermore, the extracted values are consistent with the ones obtained in Ref.[9] using 
a different method, however the present extraction method gives more accurate results. 

In the polarized case (right panel) the HTs are small within the given precision, for 
Q 2 > 1.7 GeV 2 , but they appear to drop dramatically below zero for lower Q 2 values. 
The inclusion of TMC and LxR renders these terms consistent with zero at the larger 
Q 2 values, but it does not modify substantially their behavior at lower Q 2 . From a 
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FIGURE 1. Ratio between the integrals of the measured structure functions and the calculated ones 
plotted as a function of Q 2 for the unpolarized case (left) and the polarized one (right). The calculation 
includes one by one the effects of NLO pQCD. 



comparison with results of ratio including phenomenological parametrizations 0] that 
includes some extra non perturbative behaviors it's possible to see that their effect seems 
not be large. 

As mentioned in fllOfl, an accurate extraction of the Q 2 dependence is fundamental. 
The results shown in [7] have been extendend with the recent results available in litera- 
ture for the unpolarized and polarized case. In Fig. O the higher twist coefficients of the 
present extraction in the resonance region are compared with all existing results of HT 
coefficients calculated in the DIS region. For the unpolarized case there is the BCDMS 
evaluation lfl2ll (already shown in 0]) and the new MRST calculation [11]. The HT co- 
efficients have been calculated following the factorization formula, displayed in Eq. d 
which can be expressed as fLT+ht = pLT . ( 1 +C(x) / Q 2 ). For the polarized one the only 
data available [13] are using the additive formula, for which f^t+ht _ pLT j q2 

In the expression of C(x) and H(x) there is no Q 2 dependence hidden. A different behav- 
ior for the unpolarized and polarized HT terms is evident. In details, for the unpolarized 
case in the region of high x there is a big discrepancy between the HT terms in the res- 
onance region (C res (x)) and in the DIS region (Cdjs(x)). In the polarized case, at high x 
this comparison is little bit complicated, due to the fact that there is only one point with 
x >0.6 in the resonance region (H res (x)) and no value for the DIS region (Hois(x))- 
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FIGURE 2. HT coefficients extracted in the resonance region according to Eq.([0. Shown in the figure 
is the quantity 1 +Cht(x)/Q 2 - For comparison the values for HT coefficients obtained in ref. [14] using 
DIS data and the effect of TMC are shown. The left (right) panel refers to the unpolarized (polarized) 
case. 
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FIGURE 3. HT coefficients for the unpolarized (left) and polarized (right) case. The full points rep- 
resent the values in the resonance region, while the empty points and the curves are related to the DIS 
region. 
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